Continuous-time systems

For continuous-time state equations,
continuous time state equations#Definition

Definition

Suppose homogeneous system (unforced system) ๐ฑห™=A(t)๐ฑ(t)\dot{\mathbf{x}} = A(t) \mathbf{x}(t), then, if a solution exists such that Qห™(t)=A(t)Q(t)\dot Q(t) = A(t) Q(t), then we may define transition matrix,

ฯ•(t,t0)=Q(t)Qโˆ’1(t)\phi(t,t_0) = Q(t) Q^{-1}(t)

such that

๐ฑ(t)=ฯ•(t,t0)๐ฑ(t0)\mathbf{x}(t) = \phi(t, t_0) \mathbf{x}(t_0)

Properties

  1. ฯ•(t0,t0)=I\phi(t_0,t_0) = I
  2. ฯ•(t2,t1)ฯ•(t1,t0)=ฯ•(t2,t0)\phi(t_2,t_1)\phi(t_1,t_0) = \phi(t_2,t_0)
  3. ๐ฑ(t2)=ฯ•(t2,t1)๐ฑ(t1)=ฯ•(t2,t1)ฯ•(t1,t2)๐ฑ(t2)\mathbf{x}(t_2) = \phi(t_2, t_1) \mathbf{x}(t_1) = \phi(t_2, t_1) \phi(t_1, t_2) \mathbf{x}(t_2)
    1. ฯ•(t2,t1)ฯ•(t1,t2)=I\phi(t_2,t_1)\phi(t_1,t_2) = I
    2. for all finite t1,t2t_1, t_2 on II, ฯ•(t1,t2)=ฯ•โˆ’1(t2,t1)\phi(t_1, t_2) = \phi^{-1} (t_2, t_1)

Evaluation

Given transition matrix ฯ•(t,t0)\phi(t,t_0), A(t)A(t) can be evaluated as follows

๐ฑห™(t)=ฯ•ห™(t,t0)๐ฑ(t0)\dot{\mathbf{x}}(t) = \dot{\phi}(t,t_0)\mathbf{x}(t_0)

Also,

๐ฑห™(t)=A(t)๐ฑ(t)=A(t)ฯ•(t,t0)๐ฑ(t0)\dot{\mathbf{x}}(t) = A(t)\mathbf{x}(t)= A(t) \phi(t,t_0) \mathbf{x}(t_0)

Furthermore,

ฯ•ห™(t,t0)=A(t)ฯ•(t,t0)ฯ•ห™(t,t0)|t0=t=A(t) \begin{aligned} \dot{\phi}(t,t_0) = A(t) \phi(t,t_0) \\ \dot{\phi}(t,t_0)|_{t_0 = t} = A(t) \end{aligned}

Solution of forced system equations

Consider forced system with input present,

๐ฑห™(t)=A(t)๐ฑ(t)+B(t)๐ฎ(t)\dot{\mathbf{x}}(t) = A(t) \mathbf{x}(t) + B(t)\mathbf{u}(t)

Assume complete solution has form,

๐ฑ(t)=ฯ•(t,t0)๐Ÿ(t)\mathbf{x}(t) = \phi(t,t_0) \mathbf{f}(t)

Seek conditions of ๐Ÿ(t)\mathbf{f}(t) to make this assumption true.

Thus,

๐ฑห™(t)=ฯ•ห™(t,t0)๐Ÿ(t)+ฯ•(t,t0)๐Ÿห™(t)=A(t)ฯ•(t,t0)๐Ÿ(t)+ฯ•(t,t0)๐Ÿห™(t)=A(t)๐ฑ(t)+ฯ•(t,t0)๐Ÿห™(t)\begin{aligned} \dot{\mathbf{x}}(t) & = \dot{\phi}(t,t_0)\mathbf{f}(t) + \phi(t,t_0) \dot{\mathbf{f}}(t) \\ & = A(t) \phi(t,t_0) \mathbf{f}(t) + \phi(t,t_0) \dot{\mathbf{f}}(t) \\ & = A(t) \mathbf{x}(t) + \phi(t,t_0)\dot{\mathbf{f}}(t) \end{aligned}

Therefore

ฯ•(t,t0)๐Ÿห™(t)=B(t)๐ฎ(t)\phi(t,t_0)\dot{\mathbf{f}}(t) = B(t) \mathbf{u}(t)

then

๐ฑ=ฯ•(t,t0)๐ฑ(t0)+โˆซt0tฯ•(t,ฮป)B(ฮป)๐ฎ(ฮป)dฮปy(t)=C(t)ฯ•(t,t0)๐ฑ(t0)+โˆซt0tC(t)ฯ•(t,ฮป)B(ฮป)๐ฎ(ฮป)dฮป+D(t)๐ฎ(t)\begin{aligned} \mathbf{x} & = \phi(t,t_0) \mathbf{x}(t_0) + \int_{t_0}^t \phi(t,\lambda) B(\lambda) \mathbf{u}(\lambda) d\lambda \\ y(t) & = C(t) \phi(t,t_0) \mathbf{x}(t_0) + \int_{t_0}^t C(t) \phi(t,\lambda) B(\lambda) \mathbf{u}(\lambda) d\lambda + D(t) \mathbf{u}(t) \end{aligned}

where C(t)ฯ•(t,t0)๐ฑ(t0)C(t) \phi(t,t_0) \mathbf{x}(t_0) is the zero-input response and โˆซt0tC(t)ฯ•(t,ฮป)B(ฮป)๐ฎ(ฮป)dฮป+D(t)๐ฎ(t)\int_{t_0}^t C(t) \phi(t,\lambda) B(\lambda) \mathbf{u}(\lambda) d\lambda + D(t) \mathbf{u}(t) is the zero-state response.

Fixed systems

Constant AA matrix

General solution

ฯ•(t,t0)=eA(tโˆ’t0)\phi(t,t_0) = e^{A(t-t_0)}

note, this is a matrix exponential,
matrix exponential#Definition

๐ฑ(t)=eA(tโˆ’t0)๐ฑ(t0)\mathbf{x}(t) = e^{A(t-t_0)} \mathbf{x}(t_0)

Solution of the forced system equations

๐ฑ(t)=eAt๐ฑ(0)+โˆซt0teA(tโˆ’ฯ„)Bu(ฯ„)dฯ„h(t)=CeAtB1(t)+Dฮด(t)\begin{aligned} \mathbf{x}(t) = e^{At} \mathbf{x}(0) + \int_{t_0}^t e^{A(t-\tau)}Bu(\tau) d\tau \\ h(t) = Ce^{At} B 1(t) + D \delta(t) \end{aligned}

Discrete systems

In the case of discrete time state equations,
discrete time state equations#Definition

The transition matrix takes the form

ฯ•(k,i)โ‰œA(kโˆ’1)A(kโˆ’2)...A(i)forย kโ‰ฅi+1\phi(k,i) \triangleq A(k-1) A(k-2) ... A(i) \quad \text{for } k \geq i+1

#incomplete


References

  1. https://crrl.poly.edu/6253/lectures/lect4.pdf
  2. P. E. Sarachik, Principles of Linear Systems, Cambridge Press, 1996, pp. 60-88.
  3. https://en.wikipedia.org/wiki/State-transition_matrix
  4. https://people.engr.tamu.edu/spb/courses/linearsystems/ECEN605Lecture07.pdf
  5. https://faculty.washington.edu/ratliffl/teaching/2020_LinearSystems.pdf p 12